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We construct an analytical solution to the integral equation which is believed to describe logarith- 
mic growth of the anomalous dimensions of high spin operators in planar M = 4 super Yang-Mills 
theory and use it to determine the strong coupling expansion of the cusp anomalous dimension. 
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1. Introduction: The cusp anomalous dimension is an 
important observable in four dimensional gauge theories 
ranging from QCD to maximally supersymmetric N = 4 
Yang-Mills theory (SYM) since it governs the scaling be- 
havior of various gauge invariant quantities like logarith- 
mic growth of the anomalous dimensions of high-spin 
Wilson operators, Sudakov asymptotics of elastic form 
factors, the gluon Regge trajectory, infrared singularities 
of on-shell scattering amplitudes etc. By definition 0, 0] , 
TcuspCs) measures the anomalous dimension of a Wil- 
son loop evaluated over a closed contour with a light-like 
cusp in Minkowski space-time. It is a function of the 
gauge coupling only and its expansion at weak coupling 
is known in QCD to three loops 3 and in TV = 4 SYM to 
four loops [4j. Recently, a significant progress has been 
achieved in determining r cusp (g) at strong coupling in 
planar M = 4 SYM. Within the AdS/CFT correspon- 
dence [5|, T cusp (g) at strong coupling is related to the 
semiclassical expansion of the energy of folded string ro- 
tating in AdS3 part of the target space @,0| (see also [§]) 



r cusp (5) = 2 ff - 



3 In 2 

2vr 



+ 0(1/ g) , g : 



with A = (7y M iV c being t' Hooft coupling. On the gauge 
theory side, the Bethe ansatz approach to calculating 
r c usp(ff) in the weak coupling limit was developed in 
based on integrability symmetry of planar Yang-Mills 
theory to one loop [lfj. This approach was recently ex- 
tended to all loops in planar J\f = 4 SYM. Comparing in- 
tegrable structures present on both sides of the AdS/CFT 
correspondence, an all-loop asymptotic Bethe ansatz was 
proposed in [11| . It involves a nontrivial scattering phase 
satisfying the crossing symmetry [l2} whose explicit form 
was found in [l]| . This led to an integral equation for the 
all-loop cusp anomalous dimension 1J, UM > the Beisert- 
Eden-Staudacher (BES) equation, 



K(2gt, 0) - Ag 2 / dt'K(2gt, 2gt')a(t') 



a(t) 

' ' 7u 

(2) 

with r cusp (g) = 8g 2 a(0). Here, the kernel K(t,t') 
is expressed in terms of Bessel functions, K(t, t') — 

Y.™„ i z nm (g)J n (t)J m (t')/(tt'), and its explicit form 

can be found in [14[. At weak coupling, the Neumann 
series solution to @ yields perturbative expansion of 



Tcusp(ff) in powers of g 2 which agrees with the known 
four-loop result Q. At intermediate coupling, the inte- 
gral equation ((2|) was solved numerically in [16]. The 
obtained solution for r cusp ((7) was found to be a smooth 
function of g that matches for g > 1 the string theory 
prediction (JTJ) with high accuracy. Moreover, an exact 
analytical solution to equation |J5J| in the limit g — > oo 
was recently constructed in [t7| leading to r cusp ((?) = 
2g + 0(g ), in agreement with ([1]) and with the numeri- 
cal analysis of [16( . Equation was analyzed further in 
[HI but it resisted an analytical solution so far. In paral- 
lel development, the computation of the two-loop 0(1/ g) 
corrections to the string theory prediction |T]) was initi- 
ated in [13]. Also, the result Q was reproduced [13] 
from the quantum string Bethe ansatz for a folded string 
rotating in AdSaxS 1 . 

In this Letter, we describe an approach to finding a 
strong coupling expansion of the solution to equation {2j . 
It allows us to determine exact analytical expressions for 
the coefficients in the 1/g— expansion of the cusp anoma- 
lous dimension {T]) to any desired order. 

Let us introduce two even/odd functions j±(—t) = 
±7±(t) 



-1 



t 



d(t) = 



1+ (2gt) 7 _(2 5 t) 



2gt 



2gt 



(3) 



Following [T7|], we expand j±(t) into the Bessel function 
Neumann series 



7+(*) = E fc >i(-ir +1 (2£0J 2fe (i) 7 2 fe , 

7-(<) - E fe >i(-l) fe+1 (2fc-l)J 2fc -i(t)72 fe -i 



(4) 



with the expansion coefficients ~ Jo°° dt' jt' ' Jk(t')la(t') 
(a = +/— for k = even/odd) and sign factors introduced 
for the later convenience. Substituting ([3]) into equation 
(PI) and separating even/odd in t parts, we find after some 
algebra that ([2]) is equivalent to the (infinite) system of 
equations 

7+(*) 7-(*) 



dt 

T 

dt 

T 
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e t/2g _l 
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J 2n (t)=0, (5) 



Jin— lit) — 2 < ^ n ' 1 



with n > 1. The cusp anomalous dimension can be read 
from small— t expansion 7-(t) = t r cusp (g)/(8g 2 ) + 0(t 2 ). 
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At weak coupling, one finds from ((5j> that 7_ (t) = J\ (t) - 



0(g 2 ) leading to T cusp (g) 



+ 0(g) in agreement 



with the known one-loop result [l], [|J . 

2. Exact solution: The system ([5]) has the following re- 
markable property. Introducing two even/odd functions 
r±(i) = 7±(*) T 7t W coth or equivalents 



27±(*) = 



1 — sech — 
25. 



r±(i)±tanh— r T (t), (6) 



we find from |(SJ) 

r0C eft 
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00 dt 
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r+(t) + r_(t) 
r_(t)-r+(t) 



•M*) = 0, 
Jbn— = <5ti,1 



(7) 



(with n > 1) and the cusp anomalous dimension is now 
given by r cusp ((7) = — 2<?r + (0). Here, in comparison with 
([5]), the dependence on g only resides in T±(t). 

At large g, we expect from © that the functions r±(t) 
admit expansion in the Bessel function Neumann series 



r+(t) = E fe >o(-i) fc+1 ^ fc (i)r 2fe , 



(8) 



r-(*) = J2 k >o(-i) k+x J2k-i(t)r 



2k— 1 



In distinction to the first series involves Jo(t) term 
which ensures that r+(0) 7^ 0. Indeed, for t — » one 
gets from ((4]) and ((6]) that 7+(i) ~ i 2 and 7-(t) ~ i and, 
therefore, T+(t) ~ i° and ~ t. Also, in virtue of 

J-i(t) — — Ji(t), the coefficient in front of Ji(t) is given 
by (Ti +r_i) so that it is only the sum of two coefficients 
that is uniquely defined. We make use of this ambiguity 
to choose r_i = 1. 

Substitution of © into yields an infinite system 
of finite-difference equations for the coefficients IV Ap- 
plying standard methods, we were able to construct its 
solution for T k (with k > — 1) in the following form (de- 
tailed analysis will be published elsewhere) 



r fc 



2 k 



(0) 



1 



, 9 P 
p=i ^ 



r - r (2P"l) 



(2p) 



(9) 



(v) 

where Tu are basis functions independent on g 



{p) = £(m+p-j) (p) = Hmm-j) 

2m r(m + i)r(|)' 2m ~ 1 r(m + i- P )r(i) ' 

. (10) 

and the expansion coefficients given by series in in- 
verse powers of the coupling, = J2 r >o 9~ Tc p,r- The 
sum over p in the r.h.s. of ((5]) describes the contri- 
bution of zero modes of ([7|). Their dependence on g 
is fixed by the additional condition of scaling behavior 
of 7±(f) (see Eas-lflS]) and (JTUJ) below). Knowing the 



— coefficients we can determine the cusp anomalous di- 



mension r cusp ((/) = — 2<7r_j_ (0) = 2gTo as 
00 ^ 

r cusp ( 5 ) = 2g+Y,-p-Zl 



p=l 



fr(2p- !)+ ^r( 2p -i; 
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Let us now establish the relation between the coeffi- 
cients r„ and 7„. To this end, we return to the relation 
© and apply the identities 

secht-l = E Tl > 1 (-l) n oant an , (12) 

tanhi = E„>i(-l)"a2„-it 2n - 1 , 

where a— coefficients with even/odd indices are related to 
the Euler/Bernoulli numbers, respectively. This leads to 



27±(*) = E(-!)" +1 



(13) 



n>l 



^(V2) 2 "r ± (t) T p^(t/2) 2 "- 1 r T w 



Replacing 7± (i) and F± (t) by the series dU and (HJ) , re- 
spectively, we make use of the Bessel function Neumann 
series for (t/2) p J m (t) in the r.h.s. of (|13[) to obtain 

rn m — n+1 

72™ = E E [r 2j --i^7i._ 1 + r 2j -^ >2j .], (14) 

n=l j=0 
m m—n 

72m-l = E E [^'-l K 2m-l,2j-l + r 2j ^2m-l,2j] ■ 
n=l j=Q 

Here the notation was introduced for the coefficients 

a n /g n 



K n ■ = -- 
mj 2T(n) 



(15) 



r(|(m + j + n))T(i(m - j + n)) 
T(i(m + j-n) + l)r(i(m - j - n) + 1) 



which are given by a product of two identical even/odd 
polynomials of degree (n — 1) in variables (m ± j) 



2a„/ g " 



(m + j) n +... (m—j) n + 



(16) 

with ellipses denoting terms with smaller nonnegative 
power of (m ± j). Replacing Tj in (|14|1 by their explicit 
expressions, Eqs.jl]) and (jTUJ) , we express 72m and 72m- 1 
in terms of yet unknown coefficients Cp. In particular, 
the first two terms of their large-g expansion look as 



7fc 



-A 0) + \\c-^ + c + ^ + ] +0(l/g% (17) 



with c_ = c Y + , c + = ^Cj^ + 2al and 7^ ' given i 
terms of if^, Eq.©, 7^+1 = tI°2 = h^L 7^-1 = 



3 



.,(!). 



r'" 1 and + 

i n^_ 1} ana 7 2m+ i 



.(2) 
2(r. 



72 



1 F (3) 
2 1 - 2(m 



-1) 



The 



-7 2li 

7^ 0) — term in (JTTJ) is in agreement with findings of 17 1. 

3. Quantization conditions: In our approach, the coef- 
ficients are determined from the behavior of 7 2m and 
72m-ij Eq. (|14p . at large m. To this end, we introduce 
the functions z±(x) = 72m- 1 ± 72m and examine their 
asymptotic behavior in the double-scaling limit 



to, g 



(to 



1\2 
4 



/g = fixed . 



(18) 



Employing the approach of [16| and going through nu- 
merical analysis of z± (x) , we found that in the limit (|18|) 
the solutions to @ have the following remarkable scaling 
behavior 



z+(x) 
z-(x) 





zf(x)- 


^z^(x) 




gx 


(gx)- 3/4 


zi 0) (x)- 


z {1 \x) 

f — — 




gx 



0(i/g 2 ) 
0(i/g 2 ) 



(19) 



(r) 

where the functions z± (x) (with r > 0) do not depend 
on g and have faster-than-power decrease at large x. For 
x — > 0, small— a; expansion of z±\x) runs in integer 
positive powers of x only. Indeed, matching (|17p into 
(P| we find z^O) = 1 + c+x + C>(a; 2 ) and zi 0) (a:) = 
1 + (8c_ — 3c + )x + 0(x 2 ). For x — > oo, asymptotic be- 
havior of is controlled by the coefficients . The 
quantization conditions for follow from the require- 
ment J °° dx x p z { + (x) = finite for any given p , r > 0. 

Let us start with the leading term z±(x) in the ex- 
pansion (fl"9j) . From (fl4|) and ([9]), we evaluate z±(x) = 
72m- 1 ± 72m in the scaling limit (fT8|) and find that the 
sums in (|14|) receive dominant contribution from large j. 
This allows us to substitute the K"^ ■— kernel in (fT4|) by 
its leading asymptotic behavior (flD]) and evaluate sum 
over large j in (fT4]) by integration ^ . i— > J" cZj, leading to 



7 -5/4 



2+ 



v p>0 n>l 



(20) 

where ellipses denote terms suppressed by powers of 1/g. 
Then, taking the Laplace transform w.r.t. x we obtain 



dxe x ' s z + (x) = 



fa) 



-1/4 



2g^ 

(E s % +r (f-i))(E a « s " 

\o>0 ' Si>l 



(21) 



The sum over n can be evaluated with a help of (fT2|) as 

, , v^sin(f) , 

Y,a n s n = -(1- sees + tans) = -——121.. (22) 



n>l 



sin(f + f ) 



As a function of s, it contains an infinite number of both 
poles and zeros on the real s— axis. Requiring that the 
integrals z p = / °° dx x p z + (x) should be finite for p > 0, 
we find that, firstly, J °° dx e~ x / s z+(x) is an analytical 
function of s for Re s > and, secondly, it scales at large 
s as zo — zi/s + 0(l/s 2 ). To satisfy these conditions in 
the r.h.s. of (f2Tj) . it proves sufficient to take 



5> C 

p>0 

with ct = 



T(p -i)=£- 



r(i-*) 



•0(1/5), (23) 



and £ + the normalization factor. Putting 



s = in both sides of (J23J), we get £+ = 2[r(|)] 2 . Cal- 
culating the Laplace transform dx e~ x ' s z_(x) in the 
similar manner and imposing the same conditions as for 
z+(x) we obtain the second quantization condition 



p>0 



c;r(p-|) + 2 C +( P -i)r( J3 + i)] (24) 



r 1--M 



o(i/ ff ), 



with 

(EHi, 



and 



— ) 

In comparison with 



contains the fac- 
|) that leads 



the Laplace transform of z_(a; 
tor £^° =1 a n (-s) n = V2sin(f)/sin(^ 
to (j2~4"|) . As before, putting s = in both sides of 
we fix the normalization factor £_ = j[r(i)] 2 . Then, 
expanding both sides of the quantization conditions (|23p 
and (|24|) around s = and matching the coefficients in 
front of powers of s, we determine the coefficients 
(with p > 1) to the leading order in 1/g. In this way, 

31n2 1 _ , , _ 31n2 1 
2 



+ - + 0(l/g), cl 



47T 



0(l/ff). 
(25) 

Substituting these relations into (fTTj) , we obtain r cusp (g) 
which coincides with the string theory prediction fl}. 

To calculate subleading strong coupling corrections 
to r cusp (g), or equivalently to determine the coeffi- 
cients Cp, we expand further the Laplace transforms 
/ °° dx e~ x / s z± (x) in powers of 1 / g and require each term 
of the expansion to verify the same analyticity conditions 
as the leading term. This can be done systematically by 
applying the Euler-Maclaurin formula to the sums over 
j in the r.h.s. of (|14|) . In this manner, we obtain the 
following all-order quantization conditions 



E* p 

p>0 



:o+ 



i 

H ( 

rfi- 



2tt , 



(26) 



L \4 2tt, 



52(gs)- k $(i/g), 



k>0 



p>0 



gsj p p \gs 



r i 



2tt, 



r(i 



EM" 



fe C(i/s), 



2W fe>0 
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where (1/g) = J2 r >o €k, r 9 r and the (5-independent) 
functions Qp (x) are of the form 



k, l>0 



+l Qk[ P > Q P ( X )— E 



k+lQ2l+l 
^ k,p ' 



k, l>0 



Tcusp (5 + c l) = 2.g 

- (c 5 + 23 c 2 c 3 )5" 5 



I-C2S 2 -c 33 3 - (c 4 + 2cl)g- 



c 6 



if^c 2 c 4 + 54c§ + 25c|)g- 



(C7 + C 2 C 5 + ^ C 3 C 4 + 457 C 2 C3 ) 3 - 7 



(Cg + lP C 2 C 6 



12606 _ _ 1 7916 2 , 6864 2 

■ C3C5 + -Jq- C 4 H =- C 2 C4 



29 ws^o 1 49 ^4 1 7 
„2 , Q7/1 4\ „-8 , n f„-9\ 



+4563 e 2 c| + 374 c%) + O (g~ 



where the expansion coefficients are given by 



ci 



C4 



3 In 2 

4tt ' 
21 



1 



C2 



■/3(4), c 5 



16tt 2 
43065 



K. 



C3 



C7 



2 10 7T 4 

101303055 



2 21 7T 5 



C(5), c 6 



27 
2~i%3 
1605 



2 15 7T 6 



, (28) 



C(3), 
/3(6), 



2 30 7T 7 



C(7), 



C8 



1317645 

2 22 7T 8 



/3(8), (29) 



with C(a;) the Riemann zeta function, P(x) = 
J2 n >o ( — l)™(2n + l)~ x the Dirichlet beta function and 
K = (3(2) the Catalan's constant. We verified that the 
coefficients (|29|) are in an excellent agreement with the 
numerical values obtained within the approach of [lH ] . 

The following concluding remarks are in order. 

The reason why in (|28p we expanded r cusp (<? + ci ) in- 
stead of r cusp (g) is that the c\ — dependent terms inside 
r c usp(ff) can be resummed to all orders in 1/g by simply 
replacing g — > g + c\ . This suggests to use g' = g — c\ as 
a new parameter of the strong coupling expansion. 

A distinguished feature of the series (|28|) is that the 
coefficients in front of l/g n are given by a linear combi- 
nation of the product of ((2p+l) and (3(2r) such that sum 
of their arguments equals n. Let us compare this with 
the weak coupling expansion of r cusp ((;). The latter runs 



(27) 

Explicit expressions for the coefficients Q l k follow univo- 
cally from the Euler-Maclaurin summation formula and 
they are too lengthy to present them here. For g — > oo, 
the relations (gSJ) coincide with and for €o (°) = T (g)r 
C±, Q+(0) =T(p-i) andQ-(O) = 2(p-i)r(p+i). 
Expanding both sides of (|26)l in powers of 1/g and s 
and matching the expansion coefficients, we can deter- 
mine the functions £ fc (g) and c^(g) to arbitrary order 
in 1/g. Substitution of the resulting expression for c^(g) 
into pip yields the strong coupling expansion of the cusp 
anomalous dimension. 

4. Strong coupling expansion: Solving the quantization 
conditions (|26|) . we calculated T cusp (g) to order 0(l/g 40 ). 
The first few terms of the expansion are 



in even powers of g and the expansion coefficients only 
involve products of ^—functions of both even and odd 
arguments such that the sum of their arguments equals 
the order in g 3, [il| . 

We found that, up to order 0(l/g 40 ), all expansion 
coefficients of T cusp (g) except the first one are negative. 
In addition, at large orders in 1/g, they grow factorially 
and the asymptotic expansion is not Borel summable 



k 



r(fc-|) 

(2th?) fc 



duu 1 / 2 e 
u — 2ng 



-, (30) 



with the Stieltjes integral having a pole at u = 2ng. This 
indicates that r cusp (<7) receives nonperturbative correc- 
tion ~ g 1 / 2 e~ 2jrff proportional to the residue at the pole. 

Our prediction for the cusp anomalous dimension (|28p 
relies on the strong coupling expansion of the solution 
to the BES equation ([2]). Eventual verification of |28|) 
remains a challenge for the string theory. We would 
like to mention that our result for c 2 = K/(47r) 2 is in a 
structural agreement with the (revised) two-loop super- 
string result of [l9( | and in precise agreement with a new 
superstring computation (R. Roiban and A. A. Tseytlin, 
larXiv:0709.068ll [hep-th] ) . 
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Note added: After this paper was submitted to Arxiv, 
we learned (J. Maldacena, private communication) that 
our result for nonperturbative corrections to r cusp (q) is 
in a perfect agreement with the findings of Ref. [22j. 
As was shown in (22l |. r cusp ((7), has the interpretation 
of an energy density of a certain flux configuration and, 
as such, it receives correction proportional to m 2 with 
to ~ g 1 / 4 e _7r9 being the mass gap in the 0(6) sigma 
model. 
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